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; We analyse a 2 + 1 dimensional model with charged, relativistic fermions interacting 

through a four-Fermi term. Taking advantage of its large- N renormalizability, the various 
phases of this model are studied at finite temperature and beyond the leading order in 1 /N. 
Although the vacuum expectation value (VEV) of a charged order parameter is zero at any 
non-zero temperature, the model nevertheless exhibits a rich phase structure in the strong 
coupling regime, because of the non- vanishing VEV of a neutral order parameter and due 
to the non-trivial dynamics of the vortex excitations on the plane. These are: a confined- 
vortex phase which is superconducting at low temperatures, an intermediate-temperature 
phase with deconfined vortices, and a high-temperature phase, where the neutral order 
parameter vanishes. The manifestation of superconductivity at low-temperatures and its 
disappearance above a critical temperature is explicitly shown to be due to the vortex 
confinement/deconfinement mechanism of Kosterlitz and Thouless. The ground state does 
not break parity or time-reversal symmetries and the ratio of the energy gap to T c is bigger 
than the conventional BCS value, for iV < 22. 
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I. INTRODUCTION 



Of late, there has been considerable interest in planar field theoretical models due to 
their relevance to condensed matter systems. 1 ' 2 Furthermore, the recent observation 3 of 
the renormalizability of four-Fermi couplings in 2 + 1-dimensions in large- N perturbation 
theory, 4 ' 5 has given a tremendous boost to the systematic analysis of the phase structures 
of various interacting theories. 3 ' 6 ' 7 The possibility of observing a superconducting ground 
state in some of these phases is obviously of great practical and theoretical interest. This 
is due to the fact that the conventional explanation of this phenomenon based on the 
spontaneous breaking of a U(l) symmetry may not be applicable to a lower-dimensional 
world. 

For quite some time, it has been known from the works of Hohenberg, 8 Mermin and 
Wagner 9 and Coleman 10 that breaking of a continuous symmetry (SSB) does not take 
place in 1 + 1 dimensions due to infrared divergences. This is strictly true in the ab- 
sence of long-range interactions. Similar infrared singularities also arise in 2 + 1 dimen- 
sions at finite temperature, making the vacuum expectation value (VEV) of a charged 
order parameter vanish for non-zero temperatures. It is not very difficult to identify 
the source of this problem. The low-energy field theory in the broken symmetry phase 
is governed by a massless Bose field 4> (the Goldstone mode), which in the imaginary- 
time formalism 11 can be expanded as 0(x, r) = (p n (x.)e' luJnT , where the energy takes 
on values u n = 2mr/P, (3 being the inverse temperature. The n = mode is governed 

by S[4>o] = Jq dr J d 2 xC(4>o, <90o) = P J d?xC(4>o, dfio), a purely two-dimensional action. 
Hence, it is not surprising that the afflictions of 1 + 1 dimensional theories will manifest 
themselves here. However, it is worth mentioning that a neutral order parameter (possibly 
breaking a discrete symmetry) can have a VEV in 2 + 1 dimensions at finite temperature. 3 

In light of the above constraint, to achieve planar superconductivity, one usually 
invokes a small interplanar coupling or takes recourse in unconventional mechanisms, e.g., 
anyon superconductivity. 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 19 In this paper, an alternate model 20 exhibiting 
superconductivity in the strong coupling phase due to the vortex confinement mechanism 
of Kosterlitz and Thouless (KT), 21 ' 22 is analysed in some detail via large- A?" perturbation 
theory. The model under consideration is a 2 + 1 dimensional relativistic theory with a four- 
Fermi interaction (reminiscent of the BCS theory 23 ) in the presence of a weak and external 
electromagnetic field. In the previous work 20 , the low-energy effective action up to two- 
derivatives has been explicitly computed to leading order in 1/N using two independent 
methods. 24,25 

The primary goal of the present analysis is to go beyond the leading order in 1/N 
in order to establish the renormalizability of the model and to show the absence of any 
ultraviolet or infrared divergences that could destabilize the vacuum structure. The model 
exhibits three different phases as a function of temperature and chemical potential. The 
low-temperature phase exhibits perfect diamagnetism with a pole in the current-current 
correlator because the vortices appear in tightly bound pairs by virtue of a logarithmic 
attractive potential. It is worth mentioning that on a plane the phase degree of freedom 
of a complex order parameter is multivalued, thereby giving rise to vortex configurations. 
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At a temperature T K t, the vortices are deconfined; it is explicitly shown that this is also 
the temperature where superconductivity vanishes. This result was anticipated, but not 
proven, in some recent works in the literature. 26 ' ' 28 At a still higher temperature, the 
neutral order parameter vanishes through a second order phase transition. 

The motivation for analysing the above model lies in the appearance of relativistic 
fermions as the relevant low-energy degrees of freedom in the strong coupling limit of the 
Hubbard and related models, 29 ' 30 widely believed to be of relevance to the high-T c super- 
conductors. Furthermore, four-Fermi couplings also arise naturally in the above-mentioned 
models in the presence of doping. 30 ' 31 In the absence of precise knowledge regarding the 
form and strength of these couplings, we consider a four-Fermi term of the BCS type, pre- 
serving the relativistic invariance of the original Lagrangian. The method of analysis can 
be easily generalised to other types of couplings. Finally, although not conclusive, there is 
some experimental indication as to the KT nature of the phase transition in the high-T c 
materials, 32 giving urgency to the analyses of various models exhibiting such behaviour. 

The goal of the study is to capture the generic properties of this class of models using 
the large- N renormalizability of the four-Fermi couplings in 2 + 1 dimensions and find 
the similarities and dissimilarities with the standard BCS theory. Large-iV perturbation 
theory is inherently non-perturbative: a given order in this approximation scheme repre- 
sents a class of diagrams with various powers of the conventional weak-coupling expansion 
parameter. The different orders in this perturbation theory also have nice geometrical in- 
terpretation in terms of Riemann surfaces; e.g, the leading order expression is the sum of 
all planar diagrams (genus zero surface). 5 In the present context the large- N expansion is 
ideal because of the fact that the weak-coupling perturbation is not renormalizable as can 
be easily seen from power counting. Interestingly, our analysis reveals that unlike the BCS 
case, superconductivity in the model under consideration occurs in the strong coupling 
phase. 

The organization of the paper is as follows. In Sec. II, we give a brief review of 
the derivation of the Landau- Ginzburg effective action in the BCS theory and proceed to 
do the same in our model. We start with the analysis of the effective potential V e ^ and 
show that to leading order in the 1/N expansion a neutral order parameter can have a 
non- vanishing VEV up to a critical temperature T . Subsequently, the low-energy effective 
action is computed and used to show that the VEV of a charged order parameter at any 
non-zero temperature vanishes, because of the presence of a massless mode. In Sec. Ill, 
we discuss the renormalizability of the model after taking into account the next-to-leading 
order corrections. We also show that the VEV of the neutral order parameter is not 
destabilized by quantum fluctuations. Sec. IV is devoted to the study of the occurrence 
of superconductivity in this model, taking into account the non-trivial dynamics of the 
phase degrees of freedom. Finally, we conclude in Sec. V with some comments and future 
directions of work. 
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II. ONE-LOOP EFFECTIVE ACTION 



In the discussions of BCS superconductivity by field theorists, 33 ' 34 a common approach 
is the following. One starts from the second-quantized BCS Lagrangian, including an 
effective electron-electron interaction due to phonon exchange: 

£bcs = V>} (id t - e(p)) tpi + ip\ (id t - e(p)) ^ - X^]^^. (1) 

Here e(p) is the kinetic energy, including a chemical potential and coupling to electromag- 
netism, if desired. This has the obvious U(l) phase symmetry ip — > e l6 ip. One then argues 
that the low-energy behaviour of the theory is best described in terms of a condensate field 
rather than the fermions. Physically this is connected to the fact that there is a gap in 
the fermion spectrum below the critical temperature, as a result of which long-wavelength 
excitations of the condensate are much cheaper to create than fermionic excitations. 

In order to achieve an effective theory for the condensate, one performs a Hubbard- 
Stratonovich 35 ' 36 transformation, yielding 

A = i(;\(id t - e(p))^ + ^\(id t - e(p))^ + n/WV^T + H]^\) + 9<t>*<t>- (2) 

The advantage of (2) is that its bilinear nature in the fermion fields permits one to 
eliminate the fermion, at least formally. In the language of path integrals, the generating 
functional of (2) is 

Z = J V^VijjV(j)*V(j)e lSl , (3) 

where S± = J d^x C\. Notice that the functional integral over the scalar field can be done, 
yielding an "effective action" for the fermion which is none other than the integral of (1). 
Alternatively, one can perform the fermion functional integral: 

Z = J V(j)*V(j) e iS «« 

where 



iS c ff {(t>* ,4>] 



and 



= J V^Vi)e lSl (4a) 

Here the matrix, if sandwiched between ^ = (^j,^|) and ip = (ip^ip] j ) T , gives the 
fermion-dependent part of (2). 

Next, the effective action can be evaluated in a gradient expansion. The lowest-order 
term, with no derivatives, is the effective potential for the condensate; at finite temperature, 
the result is 

Kff-alog^^^ + o^^) 2 , (5) 

1 c 
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where a, b are positive parameters, the details of which are not important to us. One sees 
from the effective potential that <fi attains a nonzero VEV below a certain critical tem- 
perature. When that happens, electromagnetism is broken and superconductivity results. 
The phase transition is of second order. 

In this paper, we pursue a similar analysis in a relativistic setting taking in to account 
the peculiarities of 2 + 1 dimensional space-time. 

Probably the first situation in which a calculation similar to the above was done in 
a relativistic setting was the work of Nambu and Jona-Lasinio. 37 This work, which is 
actually approximately concurrent with BCS, introduced the idea of dynamical symmetry 
breaking in the context of the pion-nucleon system. A four-Fermi interaction was shown 
to break chiral symmetry, leading to what are now known as Nambu- Goldstone bosons, 
which Nambu and Jona-Lasinio identified with the pions. 

Another more recent field theory which mimics the BCS theory as discussed above is 
the Gross-Neveu model, 38 a 1+1-dimensional model described by the following Lagrangian: 

£gn = tyivfcbi + y {^i^i) 2 • (6) 

Again one often introduces an auxiliary field to remove the four- Fermi interaction, replacing 
the last term in (6) with 

9<Hi^i ~ 2^ 2 - ( 7 ) 

It is well known that the discrete chiral transformation ip — > -f 5 ifj is broken spontaneously. 
This symmetry breaking, together with the fact that the model is asymptotically free, 
makes it a good testing ground for QCD (although the fact that it is in 1+1 dimensions 
makes the connection rather tenuous). 

One obvious difference between the Gross-Neveu model and the BCS theory is that 
the "condensate" <fi in (7) is a real (uncharged) scalar, while that in BCS is doubly charged. 
One can of course couple a charged scalar to fermions via an interaction of the form -i/n/>|</>| 2 , 
but this is not terribly reminiscent of BCS. In order to find a closer analog we must find 
a doubly-charged bilinear in the fermion. This can be done using the charge-conjugate 
fermion field, ifj c = Cifj T , C being the charge conjugation matrix. 

The appropriate relativistic generalization of (1), which we will analyze in detail, is 

thus 

C = ^ a (0-e4)if a --^M%^p- (8) 

Here, a, f3 range from 1 to N, N being the parameter of the large- N expansion, and re- 
peated indices imply summation. Furthermore, the coupling constant 1/A has canonical 
dimension (mass) -1 and hence conventional perturbation theory is non-renormalizable. 
However, as mentioned earlier this type of four-Fermi coupling is renormalizable in large- 
N perturbation theory, 3 making it an ideal tool to analyse the phase structure. Coupling 
to the photon field endows the model with a U(l) gauge symmetry. Like the BCS the- 
ory, it is assumed that the short range attractive force represented by the four-Fermi 
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term, originating from the underlying vibrations of the microscopic theory, dominates the 
Coulomb repulsion between the electrons at low-temperatures so that the net interaction 
is attractive. 

In what follows, we will remove the four-Fermi interaction in favour of a Yukawa-type 
coupling to an auxiliary charged scalar field <fi and then perform a derivative expansion 
of the effective theory for 0. As a first step, the effective potential will be computed to 
leading order in 1/N to study the question of radiative breaking of the U(l) symmetry d la 
Coleman- Weinberg. 39 After answering this question in the affirmative, finite temperature 
effects will then be incorporated to show that the VEV of is zero when T^O, although 
the neutral order parameter \<p\ has a non-zero VEV. At a still higher temperature the 
neutral order parameter vanishes, the phase transition being of second order. 

Were <fi a true dynamical field, rather than an auxiliary field, we would argue that the 
two-derivative term induced by the fermion loop serves only to renormalize the already- 
present kinetic term for <p. However, being an auxiliary field, it has no kinetic term and 
the induced two- derivative term is physically significant as it stands. In order for the 
expansion up to two derivatives to be sensible, therefore, the coefficient of this term must 
be positive. We check this, and find that it is. 

We use a diagrammatic approach to compute the derivative expansion, 24 making use 
of the interpretation of the effective action in terms of one-particle-irreducible Green's 
functions. The approach is essentially a generalization of the method devised by Coleman 
and Weinberg 39 to compute the effective potential. An alternative approach using an 
expansion of the trace-log representation of the effective action 25 is employed in Sec. Ill 
as a check, and also to compute the next-to-leading order effects. 

It is worth pointing out that the Lagrangian under consideration here, with additional 
scalar self-couplings, has been analysed before for the purpose of finding zero-modes of the 
Dirac equation in the presence of a vortex, either without 40 ' 41 or with 42 ' 43 a Chern-Simons 
term. Also, radiative symmetry breaking in various similar models in 2+1 dimensions has 
been discussed in Ref. 44. Hence our investigation supplements these previous works. 

Introducing auxiliary fields </> and </>* to decouple the four- Fermi term in (8), we get: 



To maintain gauge invariance of the Yukawa term in (9), the charge of the complex field 
is 2e. 

The effective action for the scalar and gauge fields is attained by integrating over the 
fermions: 



In fact, the dependence of S e g on A can be deduced from gauge invariance, so often in 
what follows we set A to zero. 



(9) 
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A nonlocal object, S e g cannot be computed exactly. However, if our goal is to describe 
the low-energy, long- wavelength excitations of the condensate, it is reasonable to attempt a 
gradient expansion of S e Q. Indeed, the equivalent of the trace-log in (4) can be computed as 
a derivative expansion directly. 25 ' 45 However, an approach which is perhaps more physical 
is to appeal to the interpretation of the effective action as an object which incorporates 
the effects of scalar field interactions mediated by internal fermion lines. Thus, S e g is the 
generating functional of one-particle irreducible vertex functions: 

S eS [0, ] = -, X; ^ J (2?r) 3 (2?r) 3 • • • (2?r) 3 (27r) 3 ( 2 ") S (5> + H> 

X r^Cpx, • • • ,p n ; q u ■ ■ ■ , q n )4>( Pl W( qi ) ■ ■ ■ 4>( Pn W(qn). (10) 

Alternatively, simply by recognizing S e s as a functional of </>, it can be expressed in terms 
of the gradients of <j) and <fi* : 

Seff[0,0*] = J d 3 x (-V^ + Zii^d^^ + Zzmrd^ + ^rf (11) 

The key feature of the first expansion is that each term is calculable: r*™'"^ is the one- 
fermion-loop diagrams with n external 0's and 0*'s. On the other hand, the second, a 
derivative expansion, is a convenient way of capturing the long-wavelength behaviour of 
the scalar field as induced by the fermion. In what follows we will compute the first 
two terms in the derivative expansion through judicious combinations of the 1PI Green's 
functions. To leading order the massless fermion propagator is z/^; the two vertices coming 
from 0, (p* couplings to fermion fields ipi and ip* are iC~ 1 5ij/2 and iC5ij/2 respectively. 

We concentrate first on the effective potential. By rewriting (11) in momentum space, 
V e ff can be calculated from the Green's functions at zero momentum since the other terms 
vanish in this limit. It is not difficult to see that 



n=l 



VM<P) = v° + tJ2 ^r^ rKn) (o, • • • 0; 0, • • • 0) . (12) 



Here, V° is the classical potential, i.e., the last term in (9). The sum over all zero- 
momentum diagrams appears formidable, but as shown by Coleman and Weinberg, it is 
actually quite doable. The only unusual twist in the present case is actually in the Feynman 
rules: the <pipip c coupling, for example, gives a vertex with an incoming <p and two outgoing 
fermions. Formally, the effective potential is given by the classical potential plus the sum 
of diagrams depicted in Fig. 1, with external fields at zero momentum. 

The above sum can be carried out, yielding 

"--^WIW 1 -^- (14) 
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V e ff is linearly divergent; we thus impose a cutoff A. Although V e g itself is not very 
illuminating, the gap equation involving vq = (\4>\), 



1 Wen v 



2Xv Q - ^ (a/A 2 + vo 2 - vq) , (15) 



N 5v Q u 2tt 

yields a non-trivial VEV for vq if A is less than a critical value A c = A/Att. Notice that 
the four- Fermi coupling is ~ 1/A and hence A < A c is a strong coupling regime. This 
result has to be contrasted with the BCS model where the gap arises in the weak coupling 
phase. The origin of this difference lies in the behavior of the density of states. In case of 
nonrelativistic fermions the density of states is taken to be constant near the Fermi surface 
whereas in the relativistic model presented here the density of states g(e) = goes to zero 
in the e — > limit. 

There are two interesting limits one can consider for the solutions of the gap equation. 
The first limit A < A c is more physical since the mean-field approach is generally sensible 
near the critical point. Here, v$ ~ 47r(A c — A); obtaining this solution we have assumed 
that vq <C A, and the solution indeed shows that this is self-consistent; consequently, the 
mass gap in the fermion spectrum is also small. Here, V e s can also be written compactly 
as 

V eS = N (A - A c ) 0*0 + ^(r<P) 3/2 , (16) 
from which it is evident that A c is indeed the critical value of A 



The second limit corresponds to the very strong coupling regime, A <C A c ; one gets 
from the solution of the gap equation vq = 2ti\ 2 /\. Here, the mass gap is large and the 
fermion mass plays the role of the cut-off. 

In what follows, we will confine ourselves to the first regime in the spirit of the mean- 
field approach. 

It is worth mentioning at this point that the presence of the linear ultraviolet diver- 
gence, although troubling at first sight, can be handled in large- N perturbation theory, 3 
where these four-Fermi couplings are renormalizable. This cannot be done in standard 
perturbation theory as seen from a simple power counting argument. In the following 
section we will give the renormalization prescription once we have the two-loop results. 

We now analyse the temperature corrections to the effective potential in order to deter- 
mine the nature of the phase transition. In the imaginary-time formalism, 11 finite temper- 
ature is achieved by making field configurations periodic or antiperiodic in time for bosons 
or fermions respectively, with period (3 = 1/T (with the Boltzman constant put to one). In 
momentum space this amounts to replacing — i J d s k/ (2n) 3 by /3 _1 J <i 2 /c/(27r) 2 , where 
k° = 2iri(3~ 1 (n + ^) for fermions. Finite density is taken into account by adding a chemical 
potential fitok . The effective potential at finite temperature and zero chemical potential 
is then 

T) = A« - i ± J |* log (l + k2 + ff + , )2 ) • (IT) 
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Like the zero temperature case, we analyse the gap equation for a non-trivial constant 
solution vt = (|0|): 

A " T "^§/(^ + 4 + W^> 2 J =a (18) 

After carrying out the sum over frequencies, we get 



- ^ .... -, 7 as) 



d 2 k 1 sinh(/? A /k 2 + ^ 

(2^)2 ^2 + ^2 cos h(/3 v ^2^T) + ! 

Introducing the cutoff A to regulate the linear divergence as for the T = case (finite tem- 
perature effects do not introduce any new ultraviolet divergence), the non-trivial solution 
v t 7^ satisfies the equation 

A - ^ + log [2(cosh(^ T ) + 1)] = 0. 

At T = 0, we get vq = 4-7r(A c — A), with A c = A/47T, which matches with our zero- 
temperature calculation. After simplification, the finite-temperature gap equation is 

ilog[2(cosh(/^ T ) + l)]=^o. (20) 



The critical temperature T is defined to be the point where = 0, yielding 

T °-2lol2- 



(21) 



In general, allowing for a finite chemical potential leads to a critical line in the \x — T plane; 
here the simplified gap equation generalizes to 

- log [2(cosh(/3u T ) + coshOfy/))] = v . (22) 

The critical line in the \x— T plane is derived from (22) by putting vt = 0, and subsequently 
solving for fx as a function of temperature. The resulting phase diagram is displayed in 
Fig. 2. 

A physically relevant quantity is the ratio of the gap parameter (at zero temperature) 
to the critical temperature. From the above, we have 2vq/Tq = 4 log 2 ~ 2.77, which is less 
than the BCS gap-to-transition temperature ratio of 3.52. However, as will be shown in 
Sec. IV, T is not the superconducting transition temperature: that temperature will be 
lower than T and hence the gap-to-transition temperature ratio will turn out to be larger 
than 2v /Tq. Indeed, it will turn out to be larger than the BCS value for < 22. 

The neutral order parameter approaches zero smoothly, with critical exponent (3 = 1/2 
(this is precisely the mean- field value), so this is a second-order phase transition. 
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We proceed to compute the kinetic term of the effective action. If we are to have a 
sensible Ginzburg-Landau theory for the condensate, the two-derivative term must have a 
positive coefficient. The computation is very similar to that of the effective potential. Here 
we present the zero-temperature case; at finite temperature, the calculation is rendered 
somewhat more complicated by the fact that Lorentz invar iance is lost; this case is discussed 
in the Appendix. 

For the kinetic term, it is clearly Green's functions with two nonvanishing external 
momenta which contribute. In real space these terms can be arranged as in (11) to exhibit 
the charge conjugation symmetry of the original theory. (This arrangement also makes 
it easier to take into account gauge interactions by minimal substitution: notice that the 
last term is already gauge covariant.) To calculate the constant coefficients Z\ and Z 2 , we 
compute the following derivatives of the two-derivative part of (11): 



5 2 S kin 



(Zi + 2Z 2 v 2 )d^5 3 



(23a) 



6(f)(x)6(f)*(0) 



<t>=4>*=v 



and 



5 2 S kin 



2Z 2 v 2 d^S 3 



(x). 



(236) 



6(f)(x)6(t)(0) 



4>=<t>*=v 



Hence, 




Z 2 d^S(x) 



1 S 2 S 



(24) 



2v 2 5(p(x)5(p(0) ' 



In momentum space, the above expressions become 




(25) 



Here, the D i *'s are the inverse propagators for the respective fields, which can be explicitly 

calculated. Consider first . It gets contributions from all one-loop diagrams where 
and cjf carry momenta p and —p, respectively; analytically, we have 




ii*{q)fc-q)D7}.{q). 
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The above series is depicted in Fig. 3. 
One finds 

1 n -i , , A -P 2 + 2v 2 J^p 2 

tt-L-'aa* (PI = ~. t; ; arctan — . 

N M KFJ 4tt 2tt 4ttv/^2 2v 

Hence, 



d 2 



dq^dqv 



q=0 



A similar calculation yields 



d 2 



dqadq 



fj, 



9=0 



5N 
167rt;o 



N 
16nv o ' 



whence 



167TW0 

AT 

Z 2 = ^7 • (26) 

1927rz;o 3 

Incorporation of electromagnetism is trivial; since we know (f> has charge 2e, the kinetic 
term must be modified to an appropriate covariant derivative; furthermore, the last term 
in (11) is already gauge-covariant. Thus, we get the zero-temperature effective Lagrangian 



£eff = 



iej~K^ + 2 * eA M 2 ~ ^* d ^ + + <pd^n - v eS . (27) 



At this point it is worth asking if a Chern-Simons term can be present in our effective 
action, thereby leading to parity and time-reversal violating effects. This term could arise 
from one-loop diagrams with two external photons carrying momentum (in addition to any 
number of scalar fields at zero momentum). However it is easy to see that such diagrams 
can never generate a Chern-Simons term. This is essentially due to the fact that there must 
be an even number of external scalar legs (to preserve gauge invariance), and therefore an 
even number of fermion propagators. Since the fermion is massless there results a trace 
of an even number of propagators; however with two-component spinors a Chern-Simons 
term would only arise from the trace of three gamma matrices (proportional to e^fc). 
Thus, there is no induced Chern-Simons term and no spontaneous parity or time-reversal 
violation, in contrast with the simplest anyon superconductivity models 46 and in accord 
with experiments on high-temperature superconductors. ' 

Once we have the full effective action, the expectation value of the charged order 
parameter can be computed at finite temperatures. Writing = pe , it will be sufficient 
to analyse the VEV of when fluctuations of p can be neglected. Hence, 

((f)(x)) = (pe id ^) ~ v T e-^ d ^ x)) . (28) 
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Going over to momentum space and expanding 6>(x, r) = J2 n 9 n (x.)e tuJnT , where u n = 
2nn/(3, we have 

<^)>4E / 7^<*»(q)*-»(-q)> 

The propagator, whose derivation can be found in the Appendix, is 

= iV„ T tanhV T /2) q2 + (i _ Z^) (30) 
Substituting this into (29), the sum can be performed, yielding 

(0 2 ) = 2 —u, r dq coth- Ul 



Nvt [l ~ ^£^) V2 tanh(^ T /2) J ° 2 (l - iT ^ y ) ^ 

The resulting integral over q = |q| diverges in both the ultraviolet and infrared; introducing 
respective cutoffs A and 77, we get 

(o 2 ) = 



Nvt (l- i 4fe y ) 1/2 tanh(^ T /2) 
up to finite parts. 

We can eliminate the above ultraviolet divergence (a next to leading order effect) by 
suitably renormalizing p, appearing in (28). To show how this happens, it is sufficient 
to work at zero temperature (since finite temperature effects do not introduce any new 
ultraviolet divergence) . The VEV of 0, then, reduces to 

A 

{4>) = v e Nv o 

which up to order 1/N becomes 

<</>> « v - A (31) 

Now recall that at the leading order vq = 47r(A c — A) with A c = After integrating out the 
9 degree of freedom, as we have done in (28) , we expect the critical coupling to get shifted, 
the shift being of order 1/N in the next-to-leading order. As will be shown explicitly in the 

following section, the new critical coupling up to order 1/N is Ac = (1 — jif)£~ an< ^ 
new VEV is (<p) = 4-7r(A c 1 ' ) — A), which precisely agrees with the expression in (31). This 
justifies our assertion that the ultraviolet divergence can be soaked up in a redifined p. 

It is not possible, however, to get rid of the infrared divergence, the dependance of 
((f)) on the infrared cutoff 77 is: 

(</>) = const ((3i]) x , (32) 

with x = 2/N((3vt tanh(/3vr/2)). We see that in the limit 77 — > 0, (</>) — > 0, since x is 
positive. This proves that 9, the massless excitation associated with the Goldstone mode, 
is responsible for the vanishing of the VEV of the charged order parameter <j> at non-zero 
temperatures. 
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III. NEXT-TO-LEADING ORDER 
CORRECTIONS AND RENORMALIZATION 

In this section, we compute the next-to-leading order corrections to ensure that the 
quantum fluctuations do not destabilize the vacuum, either due to ultraviolet or infrared 
effects. Furthermore, we show explicitly the renormalizabilty of the theory and give a 
renormalization group invariant expression for the gap. The approach used for this purpose 
is the same as the one advocated in Ref. 3, except that we have an additional massless field. 
Interestingly, the analysis of the three point vertex shows that it is not modified in the next- 
to-leading order. This is similar to the lack of radiative corrections to the electron-phonon 
vertex at low-energy, a result known as MigdaPs theorem in the literature. 49 

For calculational convenience, we use a reducible basis of 7-matrices, and a Majorana 
basis for the spinors. A common convention for the 2x2 Dirac 7-matrices is 7 = a 2 , 7 1 = 
za 3 ,7 2 = zcr 1 , where a l, s are the Pauli matrices. We can write a complex spinor in terms 
of two real ones as -0 = xi + ^X2- With = 2 + z0i and tfj c = Cifj T = Xi ~ 1X2, the 
fermionic part of the linearized Lagrangian in equation (9) (restricting for simplicity to 
one flavour for the moment) can be written as: 

C = - 0!** - z0 2 ^75^ - \N(cj)l + 2 ) - e\M75^. (33) 

Here, we have used a reducible basis involving 4x4 7-matrices 

7 M = 7m ® 

and the four-component spinor \1/ = (xi X2) 7 '• Furthermore, the 75 matrix has been defined 
as: 




In this formulation, the method of Refs. 25,45 can conveniently be used to calculate the 
derivative expansion of the effective action. For instance, using the fact that {75, 7^} = 0, 
and that 75 = 1, we obtain, assuming that the l, s are constant, 

VeS =V° + — J Tt log# + 01 - Z0275) 

= iVA0*0 + ^iV J ^3 log(p 2 - 0*0) . (34) 

This expression is equivalent to the one derived in (14) modulo a independent term, 
but appears rather easily in this approach. We have checked for consistency that the full 
covariant action of the previous section follows from this alternative method. 

To show renormalizability, the bare Lagrangian can be written (neglecting the external 
gauge field) with three substraction constants in the form: 

£ (bare) = Z 1 *0*-Z 2 0i**-zZ 2 2 *75*- XNZ 3 (0 2 + 2 ). (35) 
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Here, Z\, Z 2 and Z 3 are overall normalizations of the fields and will drop out of the 
S-matrix. 



Henceforth, we will work in Euclidean space for convenience; the Wick-rotated func- 
tional with the fermionic sources fj and r\ reads: 

Z(fj,rj) = J VVVVVfaVfoe- f d3x{c(E) -^-* v \ (36) 



where 



2 

= Z^i$ {E) ^ + Z 2 0i# + zZ 2 (/) 2 * 7 5* + XNZ 3 (J^J (<t>\ + <j>l). (37) 
Without loss of generality, we can assume that 0i takes a VEV and hence write 



01 =^(M + 1 ), 



(38) 



where M = Mq + Mi/N and Z^ = 1 + Z^/N. Since we will be computing corrections up 
to order 1/N, the parameters have been conveniently expanded in the same powers . The 
variables 4>i and ^ 2 are the fluctuations in the respective fields. 

The Feynman rules can be read off easily from the Lagrangian. Notice that to next- 
to-leading order we have the fermion mass insertion —Mi/N, given diagrammatically by 
i — | — j and the two tadpoles given by —2\NM and — 2XZ 3 (M + Mi) arising respectively 
due to the second and last term of (37). For the proof of renormalizability, it is sufficient 
to show that the full VEV of (pi and that the connected Greens functions (x)ty (y)) , 
{4>i{x)4>i{y)),i = 1, 2, are finite. 

In the momentum space, the inverse fermionic two point function is 

r«(p) = Zi^r(p), 

where 

Mi 1 

+ Mq + 

1 



i» = ^ + m + v v 



The dashed and the dotted lines in the above expression correspond to the 4>\ and </> 2 
propagators respectively. 
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The inverse mesonic two-point function is 



= N& 2 D^(p) + 2XZ 3 + 2 



+ 2 



+ 2 



(40) 



The vertex function is 



-Z 2 (l + 



N 



+ 



N 



+ finite diagrams) 



(41) 



To leading order, the gap equation (9) can be easily derived by demanding that the sum 
of the tadpoles give zero contribution. To next-to-leading order, the tadpole contributions 
diagrammatically represented as : 



gives 



Di(0) 



= -2XZ 3 M { 



1 f d 3 < 
°~2j (2^ 



(2nr°^ M) 



dD~\q-M) 



dM 



(42) 



M=M 



Here, D 1 l (q;M) and D 2 1 (q;M) are the inverse propagators for the <f>i and 2 fields 
respectively and are given by: 



and 



Dl\ q ,M) = 2X + J |^Tr 



(4 + m) (4 -ii + M) 



M M (q 2 + 4M 2 ) q 
1 1 — — ; arctan 



7r n 



(2qn) 



2M 



D^(q,M)=2X-J |^Tr 



1 



1 



M M q 
I V — arctan n 

7T 7T 27T 2M 



(4 + m) 75 (4 - 4 + m) 



75 



(43) 



(44) 
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Note that the ultraviolet behaviour of these propagators is given by: 

4 



%(9)U.^~:- (45) 



llim q^oo 

A straightforward calculation then yields: 

Mi = -2XZ 3 n - — - log(A) + finite terms . (46) 

TV 7T Z 

Now we consider the fermion self-energy correction due to the 4>i exchange; using the large 
q behaviour of the 0i propagator the divergent part can be identified: 



Taking the trace yields 



f d A q 1 
= a(0) + ia(Q)j + finite terms . (47) 



o(0)=M »/ (&?r^ log(A) ' (48) 



Similarly, after some algebra, 

5 3 g 1 p • (p — q) 



5(0) = - lim / 



p-6 J (2tt) 3 p 2 (p - q) 2 + M 2 
2 log(A). (49) 



3tt 2 



Writing the contribution of the <p2 field to the fermion self-energy as £2(2?) = 6(0) + 
i6(0)j^+ finite terms, we obtain 6(0) = -(2M /n 2 ) log(A) and 6(0) = -(2/3tt 2 ) log(A). 
The divergent part of the <pi and <p2 vertex corrections are respectively given by c(0) = 
(2/*7T 2 ) log(A) and 5(0) = — (2/n 2 ) log(A). Hence, the renormalization condition reads 

Z x - 5(0) - 6(0) = 0, 
iiM + Mi - a(0) - 6(0) = 0, 

Z 2 + c(0) + 5(0) = 0. (50) 



These conditions give 



Zi = -^log(A), 



4 
3tt 
^2 =0 

Mi =^M log(A). (51) 
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By using the gap equation for Mi, we then obtain 

2A 16M 



2AZ 3 = -— -^log(A) (52) 



7T 



and hence 

n ,r. 2A„ 1, M 16M , 
2 ^=^d-^)-^-3^1og(A). 

At this point, with all the substraction constants properly defined, the consistency of 
the renormalization procedure can be checked by computing the inverse of the connected, 
mesonic two-point function and showing that it is finite, which we now do. 

The zero-momentum inverse two-point function is 

i-°^(I) 2 ( 2a+ /p^wtW 

+ 2AZ 3 -2Mi f Tr- ' 



1 f d 3 q 



(2vr)3 (z^ + Mq) 3 
dD a dD~ 1 „ d 2 D~ l 



+ D a - 



dM dM dM 2 



(53) 



M=M 



2 J (2tt) 
here a = 1, 2. 

A rather lengthy calculation yields 

= ^^(Z 2 - Zi) + 2AZ 3 + + ^4 + finite terms . (54) 

By using the previously calculated subtraction constants Z\ , Z 2 and Z 3 the divergent part 
of the above two-point function is zero leaving only the finite parts. 

Introducing an arbitrary subtraction point fj,, we get 

z 2 =1 

2A, 1, M 16M 0l ,A, , , 

2^3=^(1- -)-^2- ^ log(-). (55) 

By demanding that the theory not depend on the scale \i, we are led to consider M = M{jj). 
In particular, since 2AZ 3 should be independent of /i, we get 

dM 16 / <9M A \ 

"^ + 3^("^r lo ^- M ) =0 ' (B6) 

yielding M(/x) ~ / u 16 / 37r27V . By substituting this expression in 2AZ 3 , we get, 

2AZ 3 _ 2 1 1 (M^A 1 - 16 ^ 

With A = 2AZ3/A and A c = ^(1 — 1/N), the /^-function can be easily computed: 
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IV. SUPERCONDUCTIVITY 



To analyse the occurrence of superconductivity in the model presented above, we study 
the low-energy effective action using a duality transformation. 28 ' 50 This is necessitated by 
the fact that, on a plane, the phase of a complex order parameter is multivalued and gives 
rise to so-called vortex excitations; a naive calculation fails to capture the contribution of 
the vortices. As will be shown below, the description in terms of the dual variables allows 
the effective separation of the single and the multivalued components of the complex field, 
and hence is quite useful in identifying the two phases: one where vortices and anti-vortices 
are tightly bound to one another, and the other where they are free. The phase transition 
between these two regimes is the celebrated Kosterlitz-Thouless transition. 21 

In this section we will show that in the low-temperature regime, where the vortices 
occur in tightly bound pairs, the relevant contribution to the path integral comes only 
from the topologically trivial piece, leading to a ground state which exhibits supercon- 
ductivity. The photon develops a mass in this phase and the current- current correlation 
function shows a massless pole, indicative of zero resistance. When the vortices are free, 
their contribution to the path integral will be seen to precisely negate the contribution of 
the topologically trivial part, destroying superconductivity. This then proves the earlier 
assertion that Tkt is indeed the superconducting critical temperature. 

We wish to analyse the long-wavelength London limit, neglecting fluctuations of \<j)\. 
We therefore write <j) = pe 10 , with p fixed at vt- As was shown in Refs. 28,50, a direct 
substitution of this change of variables into the Lagrangian (11) is inappropriate, since 
the correspondence between <j) and 6 is not one-to-one. Rather, to capture the physics of 
vortex excitations, we make the replacement <9 M 6> — > <9 M 6> — i(p*d^,(p with <p*<p = 1. The new 
variable 9 is now treated as a single- valued field and cp describes the vortex dynamics. One 
can define a vortex current by 

j; ort = (2m)-% vX d v (<p*d x <p). (59) 

This current is identically conserved, and the corresponding integrated charge gives the 
vortex number. 

The dynamics is described by the finite-temperature version of the Lagrangian (27), 
where the finite-temperature effective potential determines the value of vt (see (20)), and 
the first term of (27) determines the Lagrangian for 6>, 0. (As remarked earlier, the 
second term is independent of these variables.) The coefficient of the first term at finite 
temperature is computed in the Appendix. For static properties we can replace Z\ by 
its spatial finite-temperature counterpart Z^ given in (A. 6b), so that the Lagrangian 
becomes 

£ eff = ^ ta * h ^T/2) 4 {d , e _ ^ + 2eA , f 
It is useful to rewrite this in terms of an auxiliary current = <9 M — i(p*d^(p + 2eA fI : 

lv CeS = ~^K 3 ^ ~ ~ + 2eAn ' (60) 
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where K = vt tanh(,9w T /2)/87r. Since 9 is now single- valued, we can integrate it out, 
obtaining the constraint d^J^ = 0; this is readily solved by 

J M = ^ vX d v a\ (61) 

here a^(x) is a vector field which is defined globally. In terms of this auxiliary gauge field 
a M the new Lagrangian, after an integration by parts, becomes: 



1 



-eff 



X -fl v - 2ee^d v A x - 2na»J™\ (62) 



N AK 

where = d^a v — d v a^. Notice that the above description of J M in terms of the gauge 
field brings in an unphysical gauge degree of freedom, which must be removed by gauge 
fixing. 



is 



The integration over a M can now be performed, and the result, in the Landau gauge, 

C eS = -27T 2 NKJ™\x)^Jl m \x) - A-KeNKJ^ t9 -^e vXp d x AP - e 2 NKF^-^F^. 

(63) 

The last term is the contribution of the topologically trivial piece and is relevant for 
superconductivity. Explicitly, neglecting for the moment the second term, which describes 
the interaction of the electromagnetic field with the vortex, the current-current correlation 
function is: 

5 2 S eS /5A^x)5A u (y) = {f{x)f{y)) = -4e 2 NK - 5 3 (x - y). (64) 

This indicates a pole at zero momentum, and hence superconductivity. To see the Meissner 
effect, namely the presence of a pole in the photon propagator at nonzero monentum, we 
have to include the Maxwell term. With the addition of £ p hoton = —\F^F^ V -, the photon 
propagator does indeed have such a pole, corresponding to a photon mass of M 2 = 6 . 
Since the scalar mass near the critical point is given by M 2 = 8v 2 , the ratio M s /M v = 



e v 16 j^ can give information about the nature of superconductivity (type I or type II). 
However, without knowledge of the dimensionful coupling constants A and e (which in 
principle are determined by the underlying model), we are not able to draw any meaningful 
conclusion. 

Now we proceed to analyse the crucial question: When can one neglect the vortex 
contribution? As we will now show, at low temperatures, when the vortices are confined, 
their contribution to the F^ u d~ 2 F^ v term is zero and hence they do not affect supercon- 
ductivity. Physically, this is perfectly reasonable: assuming that the interaction of the 
vortices with the electromagnetic field is small, in the confined phase a slowly varying 
magnetic field sees pairs of tightly bound vortices with opposite vortex charges, and hence 
the net contribution from each pair vanishes, leaving untouched the contribution of the 
single- valued part. In contrast, above the KT phase transition the vortex contribution to 
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this term is crucial: it exactly cancels the contribution from the single-valued part, thereby 
destroying superconductivity. 

To see this more explicitly, consider for simplicity only the static vortex density con- 
figurations (Jj VOrt = Oand<9oJo° rt = an d let us compute the vortex contribution to the 
current-current correlation function. Knowing that the vortex current is conserved, we can 
write 

J vort (x) = p^\x) = ^m a 5(x - xa(f)) (65) 

a 

where m a is the integer-valued vorticity and x a is the position of the a th vortex. And 
making use of the green function property in two dimensions 

^ F(X) = ^ / dyF(y) log ' X - y l 
the effective action reduces to 

Sett/N = 7T K (3 ^ mam a ' log |x a - x a / 1 + 2eK(3 ^ m a / d 2 r B(x) log |x - x a | ; (66) 

a, a' a 

without the interaction term this is the action of the familiar XY model. In the static 
limit, the contribution of the vortices to the current- current correlation function is 



0'W(-q)> 



vort <5 Sefi 



<L4*(q)<L4J(-q) 
= (S ij - q 2 (p VOrt (q)p v ° rt (-q)). 



(67) 



In real space, 

P vort (q) = E m « eiq ' Xa 

a 

and hence in the confined phase 

(/ ort (q)p vort (-q)) = ]T(e^- x ^> = -2(e^ R > + (m a = 2,3. . .). (68) 

a, a' 

The expectation value of (e iq R ) can be easily computed using the action for the XY 
model, 51 yielding 

(e iq-R ) — 

\ c / q 2 ~ 2P 

and 

(P v0rt (q)/° rt (-q)) ~ (69) 

here P = nNKfi. In the confined phase P > 2, (p vort (q)p vort (-q)> ~ q p ~ 2 , wh ere p > 4. 
With (67), this implies that the zero momentum pole in (j l (q)j J '(— q)) does not get a 
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contribution from the confined vortices. This then implies that in C e s the F^d 2 F I1V 
term and hence the Meissner effect is not affected by the vortex sector. 

As an aside, it is perhaps worth pointing out that in position space 

(p vort (x)p vort (y))~|a;-y|- 2P , 

which is the well-known power law behaviour of the correlation functions in the confined 
phase. 

In contrast, at temperatures above the vortex confinement-deconfinement transition, 
when the vortices are free and the interaction between them is screened as in a plasma, we 
expect a rather substantial contribution from that sector. In this phase directly integrating 
out the variable yields: 

C=Me 2 F^F F (70) 

which is precisely of opposite sign to that of the previously computed contribution of the 
single- valued sector (c/. (63)). This substantiates the physical argument given above that 
in the deconfined phase superconductivity is destroyed and T c = Tkt- This is one of the 
main results of this paper. 

We note that the same mechanism has been invoked in other models 26 ' 27 to explain 
superconductivity. In these works the logarithmic confining potential was due to the 
exchange of massless gauge field, however the disappearance of the Meissner effect above 
Tkt has not been shown explicitly. Our analysis may also be applicable to these models; 
this is currently under investigation. 

The confining-deconfining phase transition in this model occurs at irKNj3 = 2, or, 

2 2 N v ' 

Here vt, the expectation value of \<p\ at finite temperature, is given by (20). Eqs. (20) 
and (71) can (in principle) be solved for the KT transition temperature Tkt, which as 
explained above is also the superconducting transition temperature. It is not difficult to 
see that for any iV there is a unique solution to (71); indeed, we display the left side of 
(71) as a function of inverse temperature in Fig. 4. Its intersection with the line 8/N 
gives the KT transition temperature. Clearly Tkt increases with iV; its minimum value is 
obviously fio, the inverse of the critical temperature for the transition wherein \(f>\ obtains 
an expectation value. 

We are now in a position to estimate the ratio of the gap parameter to the super- 
conducting critical temperature, 2vq/T c = 2v /Tkt, as a function of N. This ratio is a 
decreasing function of N, with minimum value 4 log 2 ~ 2.77 as N — > oo. For N small 
(more precisely, for values of N for which 8/N ^> 1), the ratio tends to the asymptotic 
form 32/iV. (Since N is integral, this asymptotic form is reasonable for N < 3.) In many 
physically interesting models 29 related to the Hubbard model, N = 2, in which case the 
above ratio is approximately 16, considerably larger than the BCS value of 3.52. In fact, 
it is found numerically that the ratio is larger than the BCS value for N < 22. 
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V. CONCLUSIONS 



To conclude, we have presented a relativistic model in 2 + 1 dimensions which exhibits 
superconductivity up to a critical temperature T c ; unlike the anyonic theories, the ground 
state does not violate parity or time-reversal symmetries. Interestingly, the gap-to-T c 
ratio is higher than the conventional BCS value. It should be emphasized that relativistic 
field theories are not uncommon in condensed matter systems; quasi-particles satisfying 
linear dispersion relations are described by such theories. Prominent examples appear 
in magnetic systems, e.g, the 0(3) sigma model descibes the low-energy modes of the 
antiferromagnetic spin system and Dirac fermions are the relevant excitations of the so 
called flux-phase ground state. These models are widely believed to be relevant for high- 
T c materials because of their connection with the Hubbard model. 

The KT mechanism of vortex binding played a crucial role in the occurrence of super- 
conductivity in our model. These vortices arise due to the multivaluedness of the phase 
degree of freedom of the order parameter. We have explicitly demonstrated how the pho- 
ton mass vanishes once the vortices unbind above some critical temperature. Concurrently, 
the pole in the current- current correlator disappears. We note that this analysis can also 
be applied to other models 26 ' 27 ' 28 ' 52 and to effective field theories describing the quantized 
Hall effect. 50 Furthermore, we remark that the relevance of the KT transitions to quasi- 
two-dimensional systems has been discussed previously; 53 such an extension to the model 
described here can also be carried out. Work along these lines is currently under progress 
and will be reported elsewhere. 

Another crucial ingredient of our analysis is the large-iV perturbation theory; this 
technique has been extensively used in condensed matter physics. In fact it has been 
recently applied to the 0(3) sigma model leading to interesting results. 54 It should be noted 
that interesting gap-to-T c ratio is reached for N < 22 in our model; it is amusing that 
multiplets of fermions appear in the low-energy dynamics of Hubbard and related models. 29 
Hence, the precise connection of the model elucidated here to that of the microscopic lattice 
Hamiltonians needs further scrutiny to establish the values of the parameters used and the 
relevant low-energy degrees of freedom. 

We acknowledge useful discussions with V. Spiridonov, D. Arovas, B. Sakita, M.B. 
Paranjape and R. Ray, and technical help from Dr. P. Tripathy in the drawings of diagrams. 
We also thank K. Shizuya for many thoughtful comments on an earliar version of this 
manuscript. R.M. thanks the Theory Group of Lawrence Berkeley Laboratory for their 
kind hospitality while this work was being completed. This work was supported in part by 
the Natural science and Engineering Research Council of Canada and the Fonds F.C.A.R. 
du Quebec. 
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APPENDIX 



Here, we present the finite temperature propagator for 9{x) when takes a non- 
vanishing VEV. Since Lorentz invariance is lost when T ^ 0, the two derivative part of 
the effective action has to be separated in to temporal and spatial parts; the coefficients 
are different at non-zero temperatures. The two derivative part can then be written as: 

£ kin = Z [ 0) d (/>d (/>*-Z[ 1) V(/>-V(/>* + Z^\(/>*d (/> + (/>d (/>^ (A.l) 

For low-energy calculations we neglect the fluctuations in the magnitude of 0, writing 
4> = (v^e 10 . Then the last two terms drop out and the kinetic term becomes 

£ kin = 4 \z[°\d 9) 2 - ^(Vfl) 2 ] . (A2) 

After performing an integration by parts: 

S(0) = J d 3 x ^9(x)D- 1 e(x) . (A3) 

Here, the inverse propagator for 6{x) is given by: 

D^(d;oj) = 2v 2 [Z^J 1 + Z<V) . (A4) 

(i) (i) 

To calculate the Z\ ;, s and Z 2 s we proceed as for the zero temperature case. We notice 
that 

^ ( ^, (0) = + Z<» V 2 - 2Z<°>„^ + 2Z<"^V 2 ) SH,) . (A5a) 

Similarly, 

S 2 S 



5</>(x)5</>(0) 



{-2Z { 2 0) v T 2 d 2 + 2Z { 2 l) v T 2 V 2 ) 5 3 (x) . (A56) 



Although we are in fact only interested in z[ l \ they are not isolated in (A. 5a) and we 



must calculate Z^ as well. 



A straightforward calculation gives 

v(0) _ at P v t ~ s inh (/3fT 
z i __iV ™ 72 



32nv T cosh 2 (/3v T /2) 
N ( , ,„ _ (3vt \ 



tanh(/^ T /2) - , L (A6a) 



16t™ t V 2cosh 2 ((3v T /2) J 

and 

Z<P = 

23 



= — — tanh(/^ T /2) . (A66) 



Similarly, 

z?> 

and 



N / (/3i> T ) 2 tanh(^ T /2) f (3v T 
192^4 V cosh 2 (^ T /2) + ' ~ _ 



2 cosh 2 (^ T /2) 



- tanh(/?ur/2) (A7a) 



Z (D = N 

2 



f3v T 



Using these expressions, we get 



2cosh 2 (/^ T /2) 



8tt 



tanh(/3w T /2) 



(A76) 



iV WT tanh(^ T /2) ( 1 _ Pvt ) U 2 + 2 



(A.8) 



Identifying the frequency u> with the Matsubara frequency u n = 2nn/ '(3, we obtain (30), 
as desired. 
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FIGURE CAPTIONS 



1. One- loop contribution to V e Q. 

2. Critical line in the (i — T plane. 

3. One-loop diagrams contributing to the two-derivative term in S e s. 

4. Graphical solution of Eq. (71). 
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